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Abstract. We have obtained a high resolution, high S/N UVES spectrum of the bright 
QSO HE 2243-6031 to analyze the damped Lya system (DLA) observed at z = 2.33. We 
measure column densities for Hi, Ni, Alii, Sin, Pn, S II, An, Cm, Fen, Mil, Znii, Cn 
Oi, and Aim; and put upper limits on the abundances of Mnii, S in, Fern, Pill, and 
Nil. The metallicity of this system is 1/12 solar at a neutral hydrogen column 
density of log iV(Hl) = 20.7. From the observed ratios [Zn/Cr] = -0.01±0.05 and 
[S/Si] = —0.06 ±0.03 we conclude that dust is very likely absent from the ISM of this 
protogalaxy. We observe an enhancement of the a/Fe-peak ratios of +0.2 dex for various 
elements, a marked odd-even effect in Mn, and a strong underabundance of N relative to Si 
and S, [N/Si,S] = — 1 at [Si/H] = —0.86. All of these ratios support an environment that is 
in an early evolutionary stage, where the onset of star formation has begun only shortly 
before the DLA was observed. We also perform a cloud-by-cloud analysis - without 
precedent at high redshift - and find a tight correlation of all low-ionization species 
with respect to Fell extending over 2.5 orders of magnitude in A(Fen). We interpret 
this trend as being due to homogeneous physical conditions (very mild ionization effects, 
common dust-destruction histories, same chemical composition) and propose that this 
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j line of sight encounters absorbing clouds that share a common environment. In addition, 

photoionization models show that these single clouds are shielded from the external 

O . ionizing radiation, so the fraction of ionized gas is small and, except for argon, does not 

Q ■ influence the measured metal abundances. The observed Alm/low-ion ratios suggest 
H ' 

' the mildly ionized gas occurs in shells surrounding neutral cores of Alii. 
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1. Introduction 

Damped Lya systems (DLAs), defined as those QSO absorption systems with the highest 
column densities, iV(Hi) > 2 x 10 20 cm~ 2 , are widely believed to occur in the ISM of the 



high redshift progenitors of present-day galaxies (Wolfe 1993). It is fortunate that such 
large column densities permit very accurate metal abundance measurements because, 
since DLAs contributed most of the neutral gas mass when the Universe was 10% of its 



present age (e.g., Storrie-Lombardi et a. 2000 ), they provide us with a powerful means 



to trace the chemical evolution of star forming galaxies over a considerable fraction of 
the Hubble time. Since the absorbing gas in DLAs is expected to be mostly neutral, the 
bulk of metals should be in neutral or singly ionized form (depending on whether or not 
the first ionization potential IP > 13.6 eV). Thus, the many transitions of these "low 
ions" that are redshifted to the optical range can be used to compute [X + /H°] ss [X/H] 
= log(X/H)-log(X/H) . 

With the advent of 8-10m-class telescopes high resolution spectroscopy of distant 
QSOs and detailed abundance studies in DLAs have become possible with relatively 
high accuracy (~ 10 %). The picture that has emerged is that of Z ~ 1/10 — 1/300 Zq 
absorbers with abundance patterns that, except for the Fe-peak ratios, resemble 



that of Milky Way metal-poor halo stars (Lu et al. 1996; Pettini et al. 1997 Prochaska 
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& Wolfe 1999; Boisse et al. 199§| ). These observations seem to support the underlying 
hypothesis that at least part of the DLAs have undergone a similar chemical evolution 
history as our own Galaxy. 

On the other hand, the observed sub-solar abundances [X/H] observed in the Galactic 
ISM are interpreted as due to a partial incorporation of the atoms from the gas phase 
into dust grains (e.g., Savage & Sembach 1996] ). Since dust is formed through conden- 
sation from the gas phase, any measure of the absorbing column density of a refractory 
element will underestimate the true number of atoms. Moreover, observations of diffuse 
clouds in the Galactic ISM have shown the clear trend of different elements to condense 
into grains in different proportions, depending on condensation temperature (Savage & 



Sembach 1996). Fe, Cr and Ni are among the most depleted elements in the ISM, while 
O, P, Ar, S, and Zn have generally very small depletion factors. Since DLAs are likely to 
sample the ISM of galaxies at earlier evolutionary stages, it seems reasonable to raise the 
question as to whether the measured abundances at high redshift have also been modi- 



fied by differential dust depletion (e.g., Lu et al. 1996 ; Pettini et al. 1997 ). Comparing 
abundances of refractory and non-refractory elements in high redshift DLAs can lead 



to estimates of the dust-to-gas ratio in those environments (Pettini et al. [19971 ; Pettini 



et al. 1999); assuming particular depletion patterns can in turn be used to correct the 



measured gas-phase abundances (Vladilo et al. 199S ) . 

However, any assessment of the relative abundances in DLAs is difficult as long as 
the interplay between dust and nucleosynthesis effects is not well understood, and may 
not be equal in high redshift DLAs and the local ISM. 

The most direct way to disentangle these two contributions to the observed abundance 
pattern is to find DLAs with as small a level of depletion as possible. Such systems, 
though, are quite scarce and only few cases are reported in the literature (Molaro et 
al. |2000| : Pettini et al. |2000| ; Prochaska & Wolfe |1999[ ). In this paper we present 
VLT UVES observations of the z = 3.01 QSO HE 2243 6031 (discovered by 



the Hamburg/ESO QSO Survey; Reimers &: Wisotzki 1997), and argue that 
the DLA observed at z = 2.33 shows negligible amounts of dust. We have 
been able to accurately measure abundances for a set of elements of different 
nucleosynthetic origins and with differing susceptibility to incorporation into 
dust grains, including the detection of the "rare" elements P and Ar. 

The paper is organized in three parts. First we address the question of the physical 
conditions in the single velocity components and perform a cloud-by-cloud analysis of 
the relative abundances (§ ||); next we look at the uncertainties introduced by possible 
ionization effects (§ ||); and the final part (§ ||) is devoted to studying dust-depletion 
effects and relative abundances. The conclusions are outlined in § 
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Table 1. UVES Observations of HE 2243-6031 



Mode Wavelength Exp. Time Observing Date 

[nm] [sec] 

Dichroic (346+580) 308-388,477-675 14 400 August 4, 5 and 6, 2000 

Dichroic (437+860) 376-500,667-1040 10 800 August 8,11 and 12, 2000 

Red arm (520) 477-676 3 600 July 10, 1999 



2. Observations and data reduction 

HE 2243—6031 was observed in service mode with the UVES instrument at the ESO 
Kueyen telescope in August 2000 under good seeing conditions (^ 0.8"). An extra 1- 
hour exposure had been previously taken during science verification of the instrument. 
Using two dichroic modes the whole optical range was covered with a total exposure time 
of 28 800 sec (Tabic 0). 

After bias-subtracting and flat-fielding of the individual CCD frames, the echelle 
orders were extracted and reduced with self-implemented routines running under MIDAS. 
The algorithm attempts to reduce the statistical noise to a minimum by fitting the 
seeing profile with Gaussians. Flux values are assigned with a variance according to 
the Poisson statistics and the read-out noise, while cosmic-ray hits are assigned with 
infinite variances. The extraction proceeds in 3 steps: (1) A Gaussian with unconstrained 
parameters is fitted to the sky profile at each wavelength using the Levenberg-Marquardt 



method (Press et al. 1986); (2) the variation of width and position with respect to the 
orders is then fitted along the dispersion direction with low-order polynomials; (3) step 
(1) is repeated, this time with width and position fixed at the values given by the fit 
solution, so only the amplitudes can vary. 

The extracted orders were wavelength calibrated using as reference Th-Ar spectra 
taken after each science exposure. All wavelength solutions were accurate to better than 
typically 1/10 pixel. The wavelength values were converted to vacuum heliocentric values 
and each order of a given instrumental configuration was binned onto a common linear 
wavelength scale of 0.043 A pixel -1 . The reduced orders were then added with a weight 
according to the inverse of the flux variances. Finally, the flux values were normalized 
by a continuum that was defined using cubic splines over featureless spectral regions. 
The spectral resolution is FWHM ~ 6.5 km s . The typical signal-to-noise ratio varies 
between S/N^ 80 at 3700-5700 A, and S/N> 50 at 5700-8000 A (becoming worse at 
both ends of the spectrum). 

In addition to the UVES data, a low resolution (FWHM « 4 A) spectrum 
of HE 2243 6031 was obtained on October 12 2001 using the Boiler and 
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Chivens (B&C) Spectrograph on the Baade 6.5m Telescope at Las Campanas 
Observatory. This spectrum was used to better define the quasar continuum 
in the spectral region around the damped Lya line. 



3. Line parameters 

Column densities N and Doppler parameters b were obtained by fitting Voigt profiles to 
each velocity component. The fits were performed using the minimum-x 2 code FITLY- 



MAN (Fontana & Ballester 1995 ). Since most of the species under study present 2 or 
more transitions, the results of the fits were in general stable and quite independent of 
the fitting region chosen. 

As an alternative approach to obtain column densities in DLAs, the apparent optical 



depth method (Savage & Sembach 1991) has been often used to get integrated column 
densities. The method is straightforward and gives accurate results for a resolved 



and non-saturated line (Lu ct al. 1996; Prochaska& Wolfe 1999). In our case, although 
most of the ions have at least one transition likely not to be saturated, we use instead the 
fitting approach, because it provides us with better column density estimates of partially 
blended components, and of those lines having evident departures from the linear part of 
the curve-of-growth. In addition, unlike previous studies, in which redshift and Doppler 
parameter of a given velocity component are held fixed for all ions, in our fitting procedure 
all fit parameters were allowed to vary.[] 

The fit results are summarized in Table ||, while the column labeled iVgt in Table || 
lists the overall column densities which resulted from summing the individual component 
column densities. All line parameters come from the Voigt profile fits, with the exception 
of iV(Znii) and iV(Crii), which were obtained using the apparent optical depth method. 
Following is a description of the fits to individual species. 



3.1. Description of the fits 

Fell. — We begin by fitting Fen because this is the ion with the most complete set of 
transitions redward of the Lya forest. We fitted simultaneously 5 Fen transitions. They 
are shown in Fig. H, plotted in velocity space with respect to z = 2.330000, and along 
with the fitted profiles (smooth line). The A2344 and A2382 transitions provided 
a fitting interval only at the position of the weaker components. We found 

1 Traditionally, line profile fitting of DLA systems has been performed by assuming 
the same Doppler parameter for all ions at given velocity component, i.e., that bulk motion 
dominates over thermal motion. This is a reasonable choice as long as DLAs occur in cold gas 
and the ions under study are massive; however, this does not prevent us a priori to chose the 
free-parameter approach described above. 
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Fig. 1. Velocity plots of observed Fen transitions (histogram) and fitted profiles. The 
dashed lines indicate the positions of the fitted velocity components. The zero- velocity 
point corresponds to z = 2.330000. 



that 14 velocity components were necessary to model the whole profiles. Their positions 
are indicated by the vertical dashed lines. We used these line positions to identify - 
but not to tie - velocity components in the remaining ions and throughout the paper 
we refer to them as to "CI", "C2", etc. We were not able to fit the absorption feature 
between Cll and C12. There appears not to be enough information on the line profile 
for F1TLYMAN to find the best line centroid. Instead of defining a fixed line position 'by 
hand' we preferred to leave the feature unfitted. An optical depth analysis shows that 
the contribution to the total Fen column density should be less than 1%. 

An. — We have detected An via the A1048 transition. The A1066 line unfortunately 
does not provide any useful constraint to the fits. Although the A1048 transition lies in 
the forest, the spectrum is clean at the expected position of its strongest components. 
Absorption by Ar I is detected in C3, C4 and Cll, which correspond to the three strongest 
lines in Fell (~ 50 % of the total Fen). The small line widths suggest these lines 
are very likely not due to weak Lya interlopers. Moreover, the position of the lines 
match well those of Fe II, with exception maybe of C4, due to blending with a Lya 
forest line (which makes the fit an overestimate of iV(Ari) in this component). This is 
the second ever measurement of argon at high redshift after the z = 3.39 detection by 
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Fig. 2. Same as in Fig. [j] but for the observed Ar I, P n and S n transitions. The dashed 
lines correspond to the Fell components. 



Molaro et al. (2001), and one of the few outside the Galaxy available (Vidal-Madjar et 



al. 2000; Levshakov, Kegel & Agafonova 2001) 



Pn. — We have detected Pn via the A963 and A1152 transitions (Fig. ||). We are quite 
confident of the detection and fit results because they are based on two transitions that 
happen to fall in a relatively unabsorbed part of the forest. This is only the third detection 



of extragalactic phosphorus so far (Molaro et al. 2001, Outram et al. 1999). The fits give 
reasonably good results for C3, C4, C6, and C8. Cll is blended with a Ly/3 line at 
z = 2.743903 with logiV = 14.57 and b = 25.1. These parameters resulted from 
fitting simultaneously the P n lines, the Ly/3, and its corresponding Lya. 



S II, Ni, Sin, Nill, Alii. — We did not encounter big difficulties by fitting these ions 
because all of them but Al n show more than one transition, although we were able to 
fit all 14 components only for Sin. The Alii A1670 line is saturated over most of the 
components, so column densities of the stronger components must be treated as lower 
limits. In Sn and Ni there is obvious blending with forest lines. CI and C2 in Ni were 
treated as one single component. 
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Fig. 3. Portion of the echelle order containing the Cm A2056 line and portion of one 
of the orders containing the Znii A2026. The flux scale was accomodated for displaying 
purposes, and the zero-velocity point corresponds to z — 2.330000. 



O I, C II. — For these ions most transitions have saturated components or they lie in the 
forest, so only a small set of components could be fitted. C3 in Cn and C8 in Oi are 
uncertain because of the small fitting interval available. CI and C2 in Oi were treated 
as one single component. 

Mnii. — Mnii is not detected (only the weak A2606 is available, the other transitions 
falling in the red cross-disperser gap at ~ 8600 A). We derive a 3cr upper limit of log N < 
11.6 for a single line. If this value is scaled according to the strongest Fell component, it 
gives a total Mnii 3<r limit of log AT(Mnii) < 12.43. 

Zn II and Cr II. — To estimate column densities for Zn n and Cr n we used the Zn n A2026 
and Cr n A2056 transitions. These are the strongest transitions available and should not 
be blended with unrelated lines. The A2062 line is partially lost within telluric lines. 

Unfortunately, the A2026 and A2056 lines lie in the noisiest part of our spectrum, 
where spurious features introduced by traces of fringing could not be removed properly. 
In addition, the Znii line happens to fall in the overlapping region of two echelle orders, 
where the continuum definition is inaccurate. This was reflected in the impossibility of 
profile fitting all components. For this reason we used the apparent optical depth method 
to constrain the respective column densities. Integrating between v = —180 and 
v = +100 km s- 1 leads to logiV(Znii) = 12.22±0.03 and logJV(Crll) = 13.24±0.02. 
We have excluded the velocity region [+100, +180], due to spurious features that 
introduce extra noise (excluding C13 and C14 introduces an uncertainty no larger than 
~ 0.5%, though). We stress, however, that due to the weakness of the lines these values 
are subject to continuum uncertainties (which are not considered in the formal errors). 

Doubly ionized species. — Fig. [7] displays the expected position in velocity space of all 
doubly ionized species covered by our spectrum. N 11 is also included in the figure. 
It is of great importance to estimate or put useful limits on column densities for these 
ions, as they will constrain the ionization level of the system. 
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Fig. 4. Same as in Fig. || but for the observed N 1 transitions. 



For Aim we have at our disposal two transitions with apparently unsaturated lines, 
so the fit results can be considered reliable. Note that C13 and C14 are not detected. 

For the rest of the ions we have attempted to put upper limits on column densities 
using different approaches. For S III we use the transition at A = 1012 A. Although the 
line lies in the forest, there is a relatively good match of the absorption profile with that 
of Fell, in line position as well as in line strengths. This suggests the S ill is present in the 
DLA gas and tracks the low-ion profiles. However, since the lines stand out only poorly 
we did not attempt a fit; instead, we used the total apparent column density, which gives 
A^(S in) = 14.05 ± 0.01 for v = [-130, +80]. 

To constrain Fein we use the A1122 transition. This line is also in the forest and 
appears so severly contaminated that no resemblance with Aim or low-ion profiles is 
recognized. An upper limit was estimated by assuming a Fein line lies at C3 and has 
b = &Fen = 8.0 km s . The maximum column density allowed by the data then becomes 
logAf = 13.4, which, weighted with C3 from the Aim fits yields a total column density 
of JV(Fem) < 14.18. Note, however, that the Fein profile does not exclude iV(Felll) = 0. 

For Pin we performed the same exercise using the expected position of the A998 
transition, which also shows no black absorption. At C3 we get logA^Pin) < 13.2 for 
b = 8.4 km s" 1 , or logiV(P in) < 13.98 cm~ 2 for the whole profile. Again, since the line is 



10 



S. Lopez et al.: A dust-free damped Lya system at z — 2.3 



12 3 4 5 6 7 8 9 10 11 12 13 14 




Fig. 5. Same as in Fig. || but for the observed O I and C n transitions. 



in the forest and we do not recognize an Al in-like absorption pattern, we cannot exclude 
AT(Pm) = 0. 

Finally, the A1083 profile, despite forest contamination, also resembles the 
low-ion pattern. This is interesting given that nitrogen is expected to be 
mainly neutral due to its high IP of 14.5 eV. Assuming that the absorption 
feature at Cll is a Nil A1083 line and weighting with Cll from the Aim fits 
we estimate that iV(Nn) < 14.55. 



3.2. Neutral hydrogen column density 

Despite its high S/N and coverage of higher order Lyman series lines, ob- 
taining AT(Hi) from the UVES spectrum alone posed some difficulties owing 
to the non-Voigt shape of the damped Lya line. To estimate AT(Hi) we used 
instead the B&C spectrum, where the fit to Lya was much better constrained 
than in the high resolution spectrum. This shows the difficulties of defining 
a continuum around an absorption feature that extends over several echelle 
orders. 

We fitted the damped Lya line with 14 individual components whose posi- 
tions were fixed at the Fen velocity components. The Hi column density was 
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Fig. 6. Same as in Fig. ||but for the observed Sill, Nill and Alii transitions. 
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Fig. 7. Same as in Fig. g but for doubly ionized species and N II. The absorption 
at v = —90 km s _1 in the top plot is unrelated from Fein 1122. 
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Fig. 8. UVES (top panels) and B&C spectrum of HE 2243-6031 at the po- 
sition of the Hi lines used to constrain iV(Hi). The smoothed line is the 
superposition of 14 Hi velocity components centered at the positions of the 
Fe II lines. The individual H I column densities are proportional to the Fe II 
column densities, and the total Hi column density log AT (Hi) = 20.67 cm 2 . 
The Doppler parameter b is assumed to be the same for all components and 
has the value of b = 16 km s \ as constrained by Ly/3 and Lye. 



varied as to keep the N(Fe n)/iV(H i) ratio constant for all components, and 
the ^-parameter was assumed the same for all components. The best fit we 
obtained implies an overall column density of logiV(Hi) = 20.67 ± 0.02 cm 2 . 
Since the same value is obtained if a single fit component is assumed, we 
conclude that the assumption of constant Fen/Hi is unimportant in so far as 
the iV(Hi) value is concerned. The fitted profiles are shown in Fig. |^. 

We note that the fit solution still leaves underabsorbed "voids" in the 
blue wings of the Ly/3 line. These are produced by low column-density Lya 
systems that do not show the low-ionization species studied in this paper. 

Of course the underlying assumption in our fitting procedure that Fen/Hi remains 
constant over all velocity components, although expected if the gas is well mixed, must 
not hold if Fe n is for instance inhomogeneously dust depleted. Next section deals with 
such possible variations. 



4. Inferences from cloud-to-cloud variations 

In this section we investigate cloud-to-cloud variations of the relative abundances. There 
are 3 possible sources of abundance variations among clouds: (1) differential depletion due 
to differing composition and destruction histories of dust grains; (2) different ionization 
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2.330613 


13.17 0.02 


5.7 0.5 




6 


2.329573 


14.01 0.00 


7.5 0.1 


Fen 


1 


2.328177 


12.51 0.016 


6.8 0.0 




8 


2.329980 


13.68 0.01 


9.1 0.3 




2 


2.328314 


12.86 0.008 


12.8 0.0 




9 


2.330252 


13.21 0.03 


10.3 1.2 




3 


2.328836 


14.09 0.003 


8.0 0.0 




11 


2.330590 


14.40 0.01 


7.0 0.1 




4 


2.329190 


14.03 0.003 


7.9 0.0 


Oi 


1+2 


2.328260 


14.18 0.01 


16.9 0.4 




5 


2.329414 


14.04 0.003 


10.6 0.0 




3 


2.328844 


15.28 0.04 


12.0 0.2 




6 


2.329595 


13.92 0.004 


6.9 0.0 




8 


2.329960 


16.15 0.30° 


7.4 1.2 




7 


2.329800 


13.60 0.005 


6.9 0.0 




13 


2.331549 


13.60 0.01 


8.9 0.4 




8 


2.329993 


13.91 0.003 


8.7 0.0 




14 


2.331820 


13.45 0.01 


6.0 0.3 




9 


2.330241 


13.45 0.005 


7.1 0.0 


Alii 


1 


2.328179 


11.76 0.03 


5.1 0.3 




10 


2.330415 


13.52 0.005 


5.7 0.0 




2 


2.328313 


12.40 0.01 


13.5 0.2 




11 


2.330608 


14.31 0.006 


4.8 0.0 




3 


2.328846 


13.03 0.01° 


8.8 0.1 




12 


2.330984 


13.20 0.005 


11.0 0.0 




4 


2.329190 


13.07 0.06° 


7.2 0.3 




13 


2.331559 


12.39 0.018 


8.2 0.5 




8 


2.329967 


13.14 0.02° 


10.5 0.5 




14 


2.331842 


12.17 0.025 


5.2 0.6 




9 


2.330226 


12.64 0.02 


6.8 0.4 


Nin 


3 


2.328828 


13.00 0.03 


5.1 0.6 




10 


2.330395 


12.56 0.02 


5.7 0.4 




4 


2.329200 


13.15 0.11 


9.5 2.6 




11 


2.330601 


13.21 0.05° 


6.3 0.3 




6 


2.329562 


13.30 0.06 


12.8 2.2 




12 


2.330989 


12.50 0.01 


15.1 0.3 




8 


2.329989 


13.04 0.10 


13.5 4.0 




13 


2.331562 


11.69 0.04 


9.4 0.8 




11 


2.330624 


13.24 0.02 


4.7 0.4 




14 


2.331798 


11.47 0.06 


11.2 1.5 


Aim 


1 


2.328165 


11.82 0.09 


7.8 2.7 


Sin 


1 


2.328164 


12.83 0.05 


4.8 0.5 




2 


2.328366 


11.82 0.08 


5.1 1.6 




2 


2.328302 


13.55 0.01 


13.8 0.3 




3 


2.328815 


12.57 0.01 


8.4 0.2 




3 


2.328842 


14.38 0.01 


8.4 0.1 




4 


2.329176 


12.47 0.01 


8.9 0.3 




4 


2.329145 


14.04 0.06 


2.9 0.8 




5 


2.329428 


12.50 0.02 


12.5 0.6 




5 


2.329301 


14.56 0.02 


13.9 0.5 




6 


2.329579 


11.67 0.11 


4.2 1.3 




6 


2.329574 


14.48 0.02 


9.6 0.4 




7 


2.329858 


12.40 0.07 


24.3 4.3 




7 


2.329815 


14.04 0.03 


7.5 0.6 




8 


2.330019 


11.68 0.09 


6.4 1.5 



14 S. Lopez et al.: A dust-free damped Lya system at z — 2.3 



1 




Sill 


i i i | i i i 
-- --x- *- 






\ 


— 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — h 
Nill 


—\ — 1 — 1 — 1 — 1 — 1 — 1 — (— 

T 


\ h- 

T 


I — I — h 


1 






— 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — h 






; 


NI 


1 1 1 | 1 1 1 1 






; 


— 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — h 
01 


H — I — I — | — I — I — I — h- 

I 


■| — h- 


I — I — b 






! 


H — l — l — l — I — l — l — l — l — I — l — l — l — l — 1 — l — l — l — l — I — l — l — l — t 
- jzn h 




i 




pii 


i i l i i i 

I 

X 












Ii " 


"f- 




bH — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — h 
SII 


H — I — I — | — I — I — I — h- 

B- 


■| — f— 

I 


I — I — E 




1 


Arl 


I I I I I I 


I 














12 12.5 13 13.5 14 

Log N(FeII) 

Fig. 9. log Ni(X)/Ni (Fell)- log Af corr (X)/iV col . r (Feil) vs. log (Fell) of all observed 
low-ion velocity components. N COTI denotes the total, corrected column densities 
listed in Table |. 

conditions due to density variations and/or local changes in the ionization field (§ |5|); 
and (3) simply differing chemical compositions. 

4-1- No abundance gradients 

Fig. |9] plots log Ni(K)/Ni(Feu) versus log TVj (Fen) for all velocity components 
for which column densities and associated errors are available. In each panel 
the values have been normalized by iV corr (X)/JV corr (FeIl), where iV corr is the total 
column density listed for each element in Table ^. Since we were not able to obtain 
individual component values of iV(Znii) we use instead iV(Feii) as reference, bearing in 
mind that iron is possibly subject to differential dust depletion. CI and C2 have been 
integrated and treated as a single component due to possible blending that could have 
affected the fits; the new value corresponds to the point at log AT (Fen) = 13.02. 
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From Fig. |9| we see that the JV(X)/iV(Feii) ratio remains, within measurement errors, 
fairly constant over almost two and half orders of magnitude in A^(Fell). This is specially 
evident in Sill, for which we have the most robust column density estimates and coverage 
of all Fe 11 components. The same trend holds even for those ions for which we were able 
to obtain N only for the strongest components (Nin, Ni, P11, S II, Arl). The major 
discrepancies with respect to the assumed abundances are found in O I and C II, but 
these are easily explained by the saturated C3 in both transitions used for C11, and by 
the overfitted C8 in O 1. Alii also shows discrepant values between weak and stronger 
components; this is most likely due to saturation effects affecting the fits to high column 
densities components (which underestimate the true AT(Aln)). 

From the data on Si 11 we deduce that the single column densities diverge from the as- 
sumed value logJV(Siil) = log A r (Fen) + 0.44 on average by only 0.15 dex. The constancy 
of AT(Siii) over such wide range of Fen column densities has important implications for 
this DLA. First, silicon is an a-element. If the delayed injection of processed material from 
SN Type I and II had undergone differing histories from cloud to cloud, one would expect 
such variations to be reflected in the individual Sin/Fen ratios; indeed, we would expect 
to observe as large a scatter in Sin/Fen as the dispersion observed among different DLAs 



of the same metallicity (e.g., Lu et al. 1996)). Secondly, silicon is a refractory element. 
Locally, the depletion of Si and Fe differ by much larger factors than what is observed 
here, over all ISM environments from the cool disk to the warm halo (Savage & Sem- 



bach 1996). In addition, since the gas-phase abundances of refractory elements decrease 
with log Hi (e.g., Wakker & Mathis fc000| ) - and thus with log AT (Fen) in this DLA- we 
would expect departures from the assumed abundances at higher values of log AT (Fell). 
Altogether, the small variations we observe in Fig. [)] mean both the chemical evolution in 
single clouds has been the same, and dust is either absent or has different properties from 
the ISM. Such small variations have been pointed out just for one other DLA (Prochaska 



& Wolfe 1996) toward Q0201+365 at z = 2.5. Those authors also suggested the absence 
of dust could explain their observations. The novelty of our approach is that we present 
the first cloud-by-cloud analysis where that trend is observed (instead of integrating the 
column densities over particular velocity windows) . From an instrumental point of view, 
the homogeneous behavior of these abundances imply the Af (Fe n)-weighted corrections 
that we will apply to some of the low-ion column densities are physically acceptable. 

The behavior we see in this DLA has not been observed in the local ISM, which shows 



a range of values of roughly 1 dex (Spitzer & Fitzpatrick 1995; Welty et al. 1999), nor 
in the LMC (Welty et al. 1999] ). It is well known that such lines of sight probe regions 
in different environments; therefore, it seems conceivable that this DLA encounters an 
homogeneous population of clouds. 
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Under the assumption that Fen tracks Hi this line of sight probes the interval 
logiV(Hi) w 18 — 20 cm" 2 , which, if in our Galaxy, would mean probing a mixture 



of halo and disk clouds (Wakker & Mathis 2000; and references therein). Wakker & 



Mathis (2000) find a remarkable anticorrelation of the metal abundance with Hi in the 
local ISM, which they interpret as the environmental dependence expected if both disk 
(with low gas-phase abundances) and halo (with high gas-phase abundances) are plotted 
together as a function of Hi. Our observations do not fit into this scenario, as - again 
assuming jV(Hi) oc TV(Fell) - they show no dependence of abundances with respect to 
Hi column density (Fig. |Io|). Two possible explanations may be responsible for this dis- 
crepancy: (1) Fe II simply does not track H I but is an increasing function of H I, in which 
case, a higher density cloud would have a higher metallicity; or (2) this particular line 
of sight probes just one type of environment. The above discussion seems to support the 
second alternative. 

5. Ionization 

If high redshift DLAs occur in dense protogalactic clumps, then one can expect a soft 
ionizing radiation field which is not able to penetrate gas with iV(Hi) > 10 20 cm~ 2 . As 
a result, all elements but those having their first ionization potential > 13.6 eV: Nl, 
Ol, An, should be primarily in singly ionized form, with little gas in ionized form (e.g., 



Viegas 1995). More recently, several authors have investigated the systematic effects of 



ionized gas on the determination of metal abundances from the low-ionization species 



e.g., Howk & Sembach 1999; Vladilo et al. 2001) 



In the present DLA the moderately low overall Af(Hi), and the detection of some 
doubly ionized ions such as Al in having similar velocity profiles as their singly ionized 
counterparts, both suggest ionization may play a role in bringing part of the low ions 
to a higher, in principle undetectable ionization level. In this section we argue against 
ionization affecting considerably our measurements of the metal abundances even for the 
single clouds. 

5.1. Self- shielding 

Our observations suggest that Hi is distributed over all of the single components. Under 
the assumption logiV(Hi) = log Af(Feii) + 5.75 for each of the velocity components, the 
median value of an Hi cloud is derived to be logiV(Hi) = 19.26, i.e., at least half of the 
clouds would be considerably ionized, were they part of a classical Lyman-limit system 



(e.g., Lopez et al. 1999). Given the complexity of the present system, however, it is also 
conceivable that self-shielding does keep the single clouds in a more neutral environment. 
Evidence for this comes from the constancy of the low-ion/Feii ratio down to the lowest 
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measurable Fell column densities. Fig. [To| shows the ratio log(iV(Siii)/iV(Feii)) as a func- 
tion of logiV(Feii). The Si n/Feii ratio of integrated column densities in DLAs should be 
only mildly influenced by ionization effects (Vladilo et al. 2001; Howk & Sembach 1999Q . 
We interpret the datapoints in Fig. |l^ as an extension of that property down to much 
lower densities, thus suggesting self-shielding in the present system. 



Next, to address the question of how ionization may influence the single-cloud abun- 
dance ratios, we compare the fit results for the single components of Fen and Aim, 
which are good proxies for the low and medium-ionization gas, respectively. As for Fe II, 
Aim column densities also fulfill the requirement of coming from unsaturated and re- 
solved lines. Ideally, we would have liked to have used Alii instead of Fell to disentangle 
the effect of possible variations in the relative abundances (despite the results shown in 
§ but the former ion has only one transition (A1670), and some of the components 
are saturated. Triply ionized species such as C iv or Si iv were not used as their profiles 
indicate they occur in distinct volumes. Fig. [ll] shows the ratio log(A''(Fell)/A ir (Allll)) as 
a function of log AT (Fen) for all 12 velocity components for which column densities and 
associated errors are available. The two lower limits correspond to the non-detection of 
Aim at the positions of C13 and C14, and the dashed line mark the Fen/ Aim calcu- 
lated from the overall column densities. We use the same y-scale as in Fig. [II] to permit 
a comparison of the scatters. Contrary to the constant behavior of the N(Si n)/N(Fe n) 
ratio, there is a larger scatter in Ar(Feii)/AT(Alm) of about 1 dex over the same inter- 
val in V(Fell). Given the constancy in the low-ion abundances, such scatter could be 
attributed to cloud-to-cloud variations in ionization. 



Our logiV(Alni)/JV(Alll) <; -0.6 is typical of DLAs with similar N(H i), according to 
the compilation by Vladilo et al. ( |2001 ). These authors find an empirical anti-correlation 
between AT (Hi) and the Aim/Alii ratio, suggesting that Hi indirectly traces the level 
of ionization in the gas. A similar trend is not seen in our Fell/Aim plot. We find 
that the two points representing the most ionized clouds (i.e., lowest Fen/ Aim ratio; 
corresponding to C1+C2 and C7) contribute only 17 % (6%) of the total Aim (Fen) 
column density. Ignoring those two datapoints we find that Fen/ Aim remains in this 
DLA quite independent of Fell (and thus of Hi in the single clouds). A linear regression 
analysis gives a slope of only m = 0.16 ± 0.32, which, should Fen indeed track Hi over 
these interval, is small when compared with the m = 0.81 ±±0.015 found by Vladilo et al. 



(2001) for larger Hi column densities. The conclusion again is that this system, although 
made of several individual absorbing clouds, can be treated as one single cloud exposed 
to the ionizing radiation. This is basically the same conclusion reached by Prochaska & 



Wolfe ( 1996 ) for the z — 2.5 DLA. Let us also note that a constancy of the ionization 
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Fig. 10. Observed ratio log N(Siu)/N (Fell) vs. logiV(Feii) along with measurement 
errors. The dashed line mark the assumed ratio from the overall Sin and Fen column 
densities. The solid curves are CLOUDY simulations for an extragalactic ionizing field 
and for [Fe/H] = -1.0 and [Si/Fe] = +0.4. 



conditions is similarly observed in the local ISM through the small variability of ionized 



helium (Wolff, Koester & Lallement |1999 ) 



5. 2. Spatial distribution of Al III 

Fig. [ll] also provides hints on the spatial distribution of Aim. In view of our finding that 
the low-ionization species in this DLA track each other independently of gas density, 
if Aim is co-spatial with the low ions then we expect that they anticorrelate, because 
a decrease in H I (low-ion) column density implies an increase in ionization. Since little 
dependence of Al in with Fe n is observed, we conclude that a scenario of partially ionized 
clouds (i.e., with Aim occuring in partially ionized interfaces bordering the neutral gas) 
is more plausible than co-spatiallity. 

In such a scenario one would expect from the line widths to observe no correlation 



between Aim and Alii 6-values. We obtain a dispersion of a = (A&) 2 /11 = 1.94 
km s _1 for Alii-Alm, while a = 1.1 km s _1 for Alu-Fen, and 0.7 km s _1 for Fell-Sill 
(including C13 and C14). If real, these differing dispersions are indicative of absorption 
occuring in different regions.n 



We stress, however, the potentially large systematic errors that may be present in our 
measurements of b, especially due to saturation of some Alii components. 
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5.3. Photoionization models 

We performed photoionization models for this DLA and compared the simulations with 
the observed Si II, Fe II, and Al ill in order to establish the level of ionization of each single 
cloud. The simulation outputs were also used to quantify possible (if any) ionization 
corrections to the observed relative abundances. 



We used the code CLOUDY (Ferland 1993) to model the ionization state for a range 



of H I column densities and ionization parameters (U, the ratio of ionizing photons to the 
total hydrogen number density). CLOUDY gives column densities for a given set of the 
input parameters: gas metallicity Z, gas number density nn, and ionizing spectrum J„. 
Of these, shape and intensity of the ionizing spectrum introduce the major uncertainties. 
We adopted the universal background ionizing field due to QSOs and galaxies calculated 
in Haardt & Madau ( 1996] ; an electronic, updated version was kindly provided by the 



authors). This is the appropriate choice for our treatment of the individual clouds as 
isolated objects embedded in an intergalactic field. Stars as sources of photoionization are 
not likely to influence significantly these models: again referring to the work by Prochaska 



& Wolfe (1996) the exponential fall off of a blackbody ionizing flux has a similar shape 
as the Hen Lyman-break in the Haardt & Madau spectrum for the energy range of 
interest. The intensity of the field at the Lyman-limit was taken to be J912 = 4.4 x 1CP 22 
ergs s _1 cm _2 Hz _1 sr _1 at z = 2.3 (consistent with results obtained using the proximity 



effect; e.g., Giallongo 1996). The geometry was plane parallel and illuminated from one 
side. The metallicity was fixed at [Fe/H] = —1.0, and the relative abundances of interest 
[Si/Fe] and [Fe/Al] were fixed at the values listed in Table ||. For fixed values of U we 
ran a grid of models by varying A^(Hi) in steps of 0.2 dex. Each integration was stopped 
when the required Hi value was reached, (i.e., the cloud size scales with Hi in order 
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to keep U constant) . The smoothed curves in figures |l0| and |ll| show the results of the 
simulations for \ogU = —2, —3, and —4. 

Fig. [l^ shows that in our model even at the very lowest ionization parameters a mild 
but significant increase of the Sin/Fen ratio with decreasing Fen is expected. Compared 
to the model, the datapoints confirm that the behavior of the observed Sin/Fen ratio is 
not consistent with isolated clouds, and that some degree of self-shielding must exist in or- 
der to keep this ratio constant over the observed range of Fe II column densities. We warn 
that the Hi value at which the iterations are stopped to obtain a given (Fen,Sin/Feii)- 
point is systematically lower than our assumed log iV(Hi) = log N (Fell) + 5.75. In other 
words, Fell is not tied to Hi in this model for AT (Fell) <^ 13.5 cm~ 2 . The simulated 
Fe II values are indeed obtained for systematically higher H I column densities than the 
assumed from the observed Fen values. However, this fact and the use in the simulations 
of slabs illuminated just from one side both support the notion that i) the individual 
clouds are embedded in a more neutral environment than if they were isolated, and ii) 
the ionization is homogeneous over the DLA components. In consequence, treating the 
system as one single cloud with N(H i) = 20.67 should not introduce major systematic 
errors. 

Concerning Al in, note that all but 2 points in Fig. |ll| lie within the contours log U = 
—4 and —3, while all Sin measurements are explained with \ogU < —4. If different 
ionization conditions explain the data, then Aim and Sin (and so Alii) may occur in 
spatially distinct regions. 

Finally, Fig. [j^ considers the same photoionization models but with AT (Hi) = 20.67. 
The vertical axis now shows the ratio N(K z )/N(X l+1 ) for the variety of observed ions as 
a function of U. The thick sections of the curves indicate those regions allowed by the 3<r 
upper limits on the doubly ionized species Aim, Fein and S III. The softer limit is set 
by Aln/Alin, but recall that Alii is probably underestimated due to mild saturation. 
Assuming \ogU « —3.8, this plot shows the approximation [X/H] w [X + /H°] 
should be accurate for the present DLA to within ~ 5 %, with the exception of 
Ar, which is still primarily ionized in this regime. 

6. Metal abundances 

6.1. Corrections to column densities 

For all ions but Fen and Sin, we were not able to determine column densities of each 
of the 14 velocity components. For weak transitions in noisy parts of the spectrum only 
the stronger components could be fitted, while for strong transitions only the weakest 
components where accessible due to blending of the saturated components. For those 
incomplete fits we applied a correction by scaling the summed column densities according 
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Table 3. Low ions: Total column densities and abundances 
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Ion 


r 




log iV c c orr 


[X/H] d 


Hi 1215 


0.41640 


20.67(0.02) 






Cn 1036 


0.12310 


>15. 02(0.06) 


>15.80 


>-1.43 0.08 


C ii 1334 


0.12780 








Ni 1199 


0.13280 


14.76(0.003) 


14.88 


-1.84 0.02 


Ni 1200.2 


0.08849 








Ni 1200.7 


0.04423 








Ni 1134.1 


0.01342 








Ni 1134.4 


0.02683 








Ni 1134.9 


0.04023 








Oi 1039 


0.00904 


16.21(0.19) 


16.79 


-0.81 0.21 


Oi 1302 


0.04887 








Oi 950 


0.00157 








Oi 948 


0.00645 








Alii 1670 


1.81000 


>13.83(0.02) 


>13.97 


>-1.18 0.04 


Sin 1808 


0.00256 


15.36(0.02) 


15.36 


-0.86 0.04 


Sin 1526 


0.12700 








Pn 1152 


0.23600 


13.28(0.08) 


13.42 


-0.82 0.10 


Pn 963 


1.46000 








Sn 1250 


0.00545 


14.88(0.01) 


15.02 


-0.92 0.03 


Sn 1253 


0.01088 








Sn 1259 


0.01624 








An 1048 


0.24400 


13.57(0.05) 


13.85 


-1.38 0.07 


Cm 2056 


0.10500 


13.24(0.02) e 




-1.11 0.09 


Mn ii 2606 


0.19270 


<11.6 


< 12.43 


<-1.77 


Fen 1608 


0.05850 


14.92(0.003) 


14.92 


-1.26 0.02 


Fen 2260 


0.00244 








Fe II 2344 


0.11420 








Fe II 2374 


0.03131 








Fe II 2382 


0.30000 








Nin 1709 


0.03240 


13.86(0.03) 


14.00 


-0.92 0.05 


Nill 1741 


0.04270 








Znn 2026 


0.48900 


12.22(0.03) e 




-1.10 0.05 



a Oscillator strengths taken from Morton (1991) for Hi, C II, Ni, Pn, SlI, An, and Mnll; 



from Verner et al. (1994) (updated values) for Oi, Alii, Aim, and S ill; from Schectman et 



al. (|1998|) for SillA1526; from Bergeson & Lawler ( |l993b| ) for SillA1808; from Bergeson & 
Lawler ( |l993a|) for Cr II and Zn II; from Fedchak et al. ( |2000|) for Ni II; from Howk et al. ( [200(]| ) 
for FellAA1608, 2260, 2344, and 2374; and from Cardelli & Savage ( |l995| ) for FellA2382. 
b Overall column densities from Voigt-profile fits. For Cn, O I, and Alii at least one compo- 
nent was saturated in all transitions available. 
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Fig. 12. Expected column density ratios as a function of ionizing parameter U for 
logiV(Hi) = 20.67 and a Haardt & Madau (|l996|) ionizing field. 



Table 4. High ions 



Ion 



r 



log iV fl t 



Aim 1862 0.28600 13.35(0.01) 6 

Aim 1854 0.57500 

S m 1012 0.04250 <14.05 d 

Fe in 1122 0.07884 <13.40 e 

P in 998 0.11200 <13.20 / 

Nil 1083 0.10310 <13.95 9 



13.35 

<14.05 
<14.20 
<14.00 
<14.55 



a Oscillator strengths taken from Morton (1991) for Fein and Nil; from Verner et al. (1994) 

for S in, Pin and Aim. 
b Overall column density from Voigt-profile fits. 

c Column densities corrected for non-fitted components using weights given by the individual 
N(Al m)-values 

d Column density limit derived from the apparent-optical-depth method 
e Column density limit derived from Voigt profile of C3 assuming bp OII = &Fem 
* Column density limit derived from Voigt profile of C3 assuming bp u = bp m 
9 Column density limit derived from Voigt profile of Cll assuming bi-ji = bun 



to the results for Fen. We have shown in the previous section that this approach is 
acceptable in this DLA. The corrected values are listed in the column labeled iV cor r of 
Table |. Since for weak transitions the strongest components dominate the total column 
densities, the corrections are small and should not introduce large uncertainties; on the 
contrary, intrinsically strong transitions with only weak components fitted (O I and C II), 
should suffer from systematic errors, and must thus be taken cautiously. 
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6. 2. Metallicity of the DLA gas 

Taking the non-refractory element Zn as reference, the gas metallicity for this DLA turns 
out to be [Zn/H] = -1.10±0.05 (or Z BhA = 1/12 Zq). A slightly lower metallicity 
is obtained if Fe is used, [Fe/H] = —1.26 ±0.02. For z = 2, this Fe abundance is quite 
at the upper end of current samples of DLA systems (Prochaska & Wolfe [2000|) , while 
the Zn value is much more representative of all measurements at z ~ 2 — 2.5 ((Pettini et 



al. |1997| ; Vladilo et al. P00l|) . This is because [Zn/Fe] > for the majority of the DLAs. 



6.3. Dust content 

The key property of the present DLA system is its apparent lack of dust. We measure 
[Zn/Cr] = —0.01 ± 0.05, a value which is consistent with no dust-depletion, since Cr, 
unlike Zn, is depleted in the local ISM. Another check for the presence of dust is provided 
by the refractory elements Fe and Ni. The Zn/Fe and Zn/Ni ratios are expected to be 
supersolar if dust is present, because Fe and Ni are two of the most heavily depleted 
elements in the ISM. In addition, [Zn/Fe] « for metal-poor stars of a wide range 



of metallicities (Gratton & Sncdcn 1991; Snedcn, Gratton & Crocker 1991: but see the 
discussion in Prochaska et al. 2000), so any departure from this value could be attributed 
to dust depletion. We observe [Zn/Fe] = +0.16 ± 0.03 and [Zn/Ni] = -0.18 ± 
0.06. While the Zn/Fe ratio implies a certain degree of dust depletion, the 
Zn/Ni abundance gives further evidence in support of dust-free gas. Altogether 
these 3 ratios indicate that the fraction of atoms that are missing from the gas phase is 
insignificant. 

Nevertheless, these inferences all rely on the goodness of our zinc abundance, which 
was not determined by profile fitting of the velocity components. An indepen- 
dent argument against significant dust depletion comes from the ratio [S/Si] 
—0.06 ± 0.03. Since (a) these two elements have a common nucleosynthetic origin, and 



(b) silicon is normally depleted in the ISM while sulfur is not (Savage & Sembach 199(i| ) , 



one expects this ratio to be supersolar if dust were present in this DLA, which is not 
observed. 

In conclusion, it seems that we have found a new DLA lacking dust, the second such 
instance at high redshift. Other cases of DLAs with small dust content are reported in 
Molaro et al. ( |2000[ ), in Pettini et al. @000| ), and in Prochaska & Wolfe (^999)) . From the 



Pettini et al. compilation the DLA with lowest dust content has [Zn/Cr] — —0.13 ±0.20 
at z = 0.86 toward Q0454±039, and from Prochaska & Wolfe, [Zn/Cr] = -0.20 ±0.12 in 
the z = 1.920 DLA toward Q2206-199. In both cases, there appears to be some degree of 
dust depletion, at least stronger than in the present DLA. A more similar case as ours 



was reported by Molaro et al. (2000) who measured [Zn/Cr] = —0.06 ± 0.08 at z = 3.39 
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toward Q0000-2620. This shows that the present system is uncommon among 
the few documented low-dust DLAs. 



6.4- a/Fe-peak ratios 

A diagnostic tool to probe the evolutionary stage in DLAs is the abundance ratio of a 
to Fe-peak elements. Due to their different origins, a-capture elements being produced 
mainly in SNe Type II, and Fe-peak elements being delivered to the ISM by SNe Type la 
on a longer timescale, the a to Fe-peak ratio is expected to rapidly attain supersolar values 
in a chemically less mature system. Moreover, this ratio should decrease as cosmic time 



proceeds, unless possible environmental factors might influence it (Ellison & Lopez 2001). 
In our Galaxy the values range from < [a/Fe] > ss +0.5 in metal-poor stars down to 
k 0, as metallicity increases. 

All a elements we measure in the DLA toward HE 2243—6031 follow each other quite 
closely (with the exception of Ar; see below). We obtain [Si/Zn] = +0.24 ±0.05, and 
[S/Zn] = +0.18 ±0.04 (recall that silicon is normally depleted in the local ISM; Savage 



& Sembach 1996 ). Even the a-element oxygen, despite the larger measurement errors, 
fits well into the trend with [O/Zn] = +0.29 ± 0.22. On average, we obtain for Si 
and S < [a/Zn] >= +0.21 ± 0.06 (simple mean). This value is in good agreement with 



measurements of metal-poor halo stars (e.g., Ryan et al. 199E and references therein) or 



the thick disk (Prochaska et al. 2000) at [Fe/H] = —1.2, but is high when compared with 
the current sample of DLAs. The conclusion is: like in Galactic metal-poor stars, these 
ratios probably reflect the temporal delay between SNe of Type I and II. They show that, 
despite its high metallicity, this DLA gas has not undergone significant metal pollution 
from SN Type I yields. We re-emphasize that these ratios are likely unbiased from 



dust-depletion effects. Centurion et al. (200C) have used the [S/Zn] ratio to investigate a 
possible evolution of the a/Fe-peak abundance ratio that is unbiased from dust effects. 
Although their sample is small (6 measurements) those authors do find a decrease in 
that ratio with increasing metallicity, which is what one would expect if metallicity keeps 
track of chemical evolution in DLAs and if the objects conform a chemically homogeneous 
sample. Our measure of [S/Zn] does not fit in that trend, being too high for [Zn/H] 
= —1.1, but it much better conforms to the enhancement observed in metal-poor stars. 



6.5. The odd-even effect 

The odd-even effect, that is, the underabundance of odd-i? elements relative to even-i? 
elements of the same nucleosynthetic origin, is another known property of Halo stars. 
Since Fe is more prone to dust depletion than Mn in a variety of ISM environments, 
this ratio can be used as a discriminant between dust depletion and pure SN Type II 
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enrichment. In fact, [Mn/Fe] averages « —0.3 in Galactic Halo stars at [Fe/H] = 
(Ryan et al. |1996[ ), whilst [Mn/Fe] « +0.3 in the local ISM (Savage & Sembach p.996|) . 

The Fe il-corrected value of our 3cr upper limit for the non-detection of Mn n yields 
[Mn/Fe] < —0.51, which confirms the odd-even effect, is in accordance with the sample of 



Halo metal-poor stars in Ryan et al. ( 1996 ), but is quite low for that of thick disk stars in 
Prochaska et al. (200C). The low value also conforms with previous DLA measurements 



(Lu et al. p_996[ ). 

We do not observe the same effect for P, the other odd-element observed in this 
DLA. We obtain [P/Si] = +0.10 ± 0.09 at [P/H] = -0.8. This value is high if 
compared with [P/Si] = -0.40 at [P/H] = -2.3 (Molaro et al. |200l|) and [P/Si] 



-0.30 + 0.09 at [P/H] = -1.2 (Outram et al. 1999). 



6. 6. Under abundant Ar I 



Our column density estimate for An implies a significant underabundance of Ar with 
respect to other a-chain elements. This corresponds to ~ —0.5 dex of what is expected 
if argon is overabundant by +0.2 dex (assuming argon tracks a elements [Timmes, 



Woosley & Weaver 1995 ). In the local ISM, An is rarely found or is significantly below 
its cosmic abundance relative to Hi (Sofia & Jenkins 1998] ; Jenkins et al. 200C). As 
suggested by Sofia & Jenkins in order to explain their finding of a low An abundance 
in the local ISM, the underabundance can be explained if a significant part of the argon 
is ionized. In this DLA that fraction is ~ 60 %. Moreover, since An has a much larger 
ionization cross section than Hi, our measurement can also be explained in terms of 
moderately ionized gas, which brings An into Aril in higher proportions than for the 
rest of the low-ions. The alternative explanation, a low An abundance caused by dust 
depletion, is in our DLA ruled out given (1) the scarce propensity of Ar to dust-grain 



incorporation (Meyer et al. 1991; Sofia & Jenkins 1998) and (2) the negligible amount of 



dust (6.3). Thus, these conditions resemble very much those in the local ISM described by 
Sofia & Jenkins where argon is systematically found underabundant, and confirm their 
suggestion that ionization, and not dust, is the responsible agent. 



6.7. Nitrogen 

Chemical evolution models predict the production of nitrogen in the CNO cycle in stellar 
interiors to have two different origins. While primary N is produced in intermediate-mass 
stars (3-8 Mq) without dependence on the initial metal content, secondary nitrogen 
is also produced in stars of any mass but with an initial supply of heavy elements. The 
net result is that the N/O ratio should be independent of O/H for primary N while 
N/O oc O/H for secondary N. In addition, since intermediate-mass stars evolve on much 
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larger scales than the short-lived progenitors of SN Type II - responsible for the bulk 
of a-elements like O, Si, S, etc. -, there should be a delay of ~ 5 x 10 8 yr between the 



release of primary N and that of O (Lu et al. 1998 and references therein). 

We measure [N/O] = —1.03 at [O/H] = —0.81. However, despite the accurate N 
abundance we were able to obtain, our oxygen estimate is rather uncertain.^ Relying 
rather on the better estimates for Si and S, we find [N/Si] = —0.98 at [Si/H] = —0.86 
(and basically the same values are found if S is used). Referring to Fig. 1 in Lu 



et al. ( |1998|) , this value is at the very low end of DLA measurements at this Si abun- 
dance; moreover, according to the fits to Galactic measurements, it is better explained 
as secondary nitrogen^. This is in agreement with the high a/Fe-peak ratios we observe 
because both effects can be explained by an early phase in this protogalaxy's chemical 
evolution, observed at a time when the release of Fe-peak elements and primary nitrogen 
has not yet begun. 



7. Summary of the conclusions 

1. We have studied in detail a z = 2.33 DLA observed at high resolution toward the 
bright QSO HE 2243—6031. We have been able to detect and obtain reliable column 
densities for H I, N I, Si II, P II, S II, Ar I, Cr II, Fe II, Ni II, Zn II, and Al III; less accurate 
determinations for Cn, Alu and Oi; and to put upper limits on Mnii, S in, Fein, 
Pin, and Nil. The metallicity of this system is [Zn/H] = -1.1. 

2. By comparing column densities of single velocity components we have performed the 
first cloud-by-cloud analysis in a high-z DLA so far. We find a tight correlation of all 
low-ionization species with respect to Fen extending over 2.5 orders of magnitude in 
A^(Fen). We interpret this trend as a lack of significant changes in the dust, ioniza- 
tion, and intrinsic chemical-enrichment effects. If Fen tracks Hi, such homogeneity 
differs from what is observed in the local ISM, meaning that our line of sight probably 
samples clouds located in one single environment. 

3. We find a probable scenario for the occurrence of Aim is ionized shells surrounding 
neutral cores of Alii. Photoionization models show that the ratio of Aim to low ions 
can be modeled with a wider range of higher ionization parameters than the ratio 

3 Fig. [] shows C8 departures from the trend followed by the other 4 components. This might 
be due to an overestimate of C8 (which dominates the total column density). Excluding C8 
and re-scaling, however, implies lowering the total column density by 0.7 dex. Although such 
value conforms better with the rest of the components, it translates into [O/H] = —1.5, so the 

underabundance of N over O would become only [N/O] = —0.2. 

4 But let us note that N II could contribute with ~ 0.20 dex to the nitrogen abun- 
dance, thus bringing [N/Si] more in line with the Lu et al. values. 
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among low ions. Independently of this, the disparate line width differences between 
Alii and Aim also support the notion of separate Alii and Aim regions. 

4. Using different abundance ratios between ISM refractory and non-refractory elements 
we have shown that the observations are compatible with negligible amounts of dust in 
this protogalaxy, although the measured [Zn/Fe] = +0.16 ±0.03 may indicate 
a small degree of dust depletion. 

5. We observe - unbiased by dust depletion effects - an enhancement of the a/Fe-peak 
ratios of +0.2 dex for various elements, a marked odd-even effect in Mn, and a 
strong underabundance of N relative to Si and S, [N/Si,S] = —1 at [Si/H] = —0.86. 
All of these ratios support an environment that is in an early evolutionary stage, 
where the onset star formation has begun shortly before the DLA was observed. The 
absence of dust is in agreement with this notion, if there has not been time enough 
to form dust grains. 

6. The observed underabundance of An by 0.5 dex is easily explained by a significant 
fraction of argon residing in a higher ionization phase. This situation resembles the 
conditions in the local ISM. Our ionization simulations demonstrate, however, such 
an ionized phase does not contain a considerable part of the other elements, so our 
abundance measurements are not affected by undetected ionized atoms. 



7.1. Final comment 

The observed uniformity of metal abundances across the line of sight to HE 2243—6031 
has important implications for current models of galaxy formation. As mentioned before, 
it suggests the individual components are embedded in the same environment, which, in 
addition, must have undergone significant mixing. If DLAs are formed from the accre- 



tion of protogalactic small clumps (e.g., Haehnelt et al. pL998| ) that have formed indepen- 
dently, then one expects each protogalaxy to contribute with gas of different dust and 
nucleosynthetic compositions. The homogeneity of metal abundances observed in this 
DLA is difficult to explain in this scenario, unless the mixing of metals has occurred in 
short enough timescales to counteract the enhancement of any original differences (e.g., 



through the metal release by SN Type I; see the discussion in Prochaska & Wolfe 2002] ) 



The chemically young gas we apparently observe might be a consequence of such a rapid 
process. 

Although the statistical weight of our result on single-cloud abundances may not 
be sufficient to draw any general conclusions on the physical nature of DLAs, the ob- 
served uniformity in HE 2243—6031 seems to be a quite common property among DLAs 
(Prochaska et al. |200l|) . 
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On larger physical scales, also important for constraining the models is to investigate 
abundance variations among systems. There are few known 'DLA groups' along single 
lines of sight where the redshift separation is small enough to support a physical connec- 
tion between the DLAs, but big enough to characterize the systems as separate entities 



(i.e., a few 1000 km s ; Lopez et al. 2001; Ellison & Lopez 2001). In the only such 



groups so far investigated at high resolution, Ellison & Lopez (2001) find very similar- 
relative abundances among the DLA members, despite their different metallicities. While 
it is not clear whether such DLA groups later merge, environment may be an effective 
catalyzer of metal enrichment in DLAs at early evolutionary stages. 
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